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Abstract 

c/j We have developed a set of college level, table-top labs that can be performed with an interfer- 

O 

ometer using satellite TV electronics and compact fluorescent lamps as microwave signal sources. 
This interferometer, which was originally developed at the MIT Haystack Observatory as a Very 

ft 

Small Radio Telescope (VSRT) to observe the Sun, provides students with hands-on experience 

in the fundamentals of radio interferometry. These labs are easily performed and convey an intu- 

y— I itive sense of how combining the signals from an array of antennas reveals information about the 

00 

structure of a radio source. 

ON 

We have also developed a package of java programs, called "VSRTLPlotter" , which is available 
H as a free-download, to facilitate the data processing and analysis of these labs. 

X 
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I. INTRODUCTION 



The radio astronomical technique of aperture synthesis has been used for extensive as- 
tronomical studies for decades. Yet, because of the complexity of the math involved, this 
important technique is almost always excluded from an undergraduate curriculum in physics 
and astronomy, leaving the interested students mystified about how instruments like the Very 
Large Array in New Mexico^ actually work. 

We present here a discussion of a set of labs for the undergraduate level which provide 
first-hand experience with the basics of making images with aperture synthesis. Although 
a thorough mathematical explanation of the process would be helpful for the advanced 
students, these labs are designed to be stand-alone so that they can be completed without 
the mathematical lectures and still impart a general conceptual understanding to the lower 
level students. 

II. SOME BASICS OF APERTURE SYNTHESIS 

A number of antennas arranged in an array with assorted baselines receive the radiation 
from a given source simultaneously, and the output signals from each pair of antennas 
are cross-correlated. Correcting for a few systematics, the cross-correlations lead to what 
is known as the "Visibility Function," which is a complex-valued function of the baseline 
vector. An image of the source is then obtained via the Fourier transform of the Visibility 
function. 

The Visibility for any given baseline vector between antennas contains an amplitude and 
a phase and so is simply represented by 

V(b) = Ae*, (1) 

where b is the projection of the baseline onto the sky plane. For the purposes of this paper, 
we can simplify this discussion and consider just the two-dimensional situation in which the 
source structure and the antennas are located in a single plane, in which case the baselines 
and the source structure are each one-dimensional. 

For a number of point sources, the visibilties due to all the individual sources simply add, 
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so that the total Visibility for baseline b is 



V T (b) = Ae l(t> = ^ A k e^ k 



(2) 



k 



The magnitude of any particular Visibility, given as the square root of the sum of the squares 



Since the Visibility function is the Fourier transform of the image, the shape of the 
Visibility function contains information about the structural aspects of the source. For 
example, if observing a pair of equally bright point sources, the Visibility function has nulls 
whose spacing in projected baseline is inversely proportional to the angular distance between 
the point sources. More specifically, the nulls are at 



where N is an odd integer, b is the projected baseline and A9 is the angular separation of 
the point sources. 

And, if observing a single but resolved source, the Visibility function will be a decreasing 
function with baseline length and the rate of decrease will be inversely proportional to the 
angular size of the source. For a source with a Gaussian brightness profile, for example, the 
Visibility function will also decrease with a Gaussian profile, and the half-maximum width 
of the Visibility function will be inversely proportional to the half-maximum width of the 
source brightness distribution. 

In the labs we discuss here, the students obtain a first hand exposure to the nature of 
the Visibility function and how it relates to the source structure. Without any discussion of 
Fourier transforms, or even complex numbers, students gain an intuitive understanding of 
how an array of radio antennas produces images of radio sources. 

III. THE VSRT INTERFEROMETER 

The labs we discuss here use the Haystack VSRT (Very Small Radio Telescope) 
Interferometer. This instrument costs less than $500, is easy to assemble, is sta- 
ble and reliable, is easy to operate, and is easily manipulated in the lab room. In- 



of the total real part and the total imaginary part, then, is given by^ 




(3) 




(4) 



3 



formation about purchasing and assembling the VRST interferometer can be found at 
http: / / www.haystack.mit .edu / edu / undergrad / VSRT / index.html 

The VSRTs come with satellite TV dishes, which are needed for observations of the Sun, 
as described by Doherty, Fish, and Needles (2011) 3 . But, for the labs in the classroom the 
dishes are not needed and so in all the labs we discuss here the feeds act as the antennas. 
In the following, then, the words "feeds" and "antennas" are interchangeable. 

The radio sources to be used in the lab room with the VSRT interferometer are, simply, 
compact flourescent light bulbs (CFLs). In addition to visible light, these light bulbs emit 
microwaves, which can be detected by the VSRT feeds. This radio emission is due to 
bremsstrahlung radiation emitted by hot, free electrons in the plasma, produced by the gas 
discharge, when they collide with the glass walls. 



A. How the VSRT Interferometer Works 



The VSRT interferometer differs from modern interferometers used by radio astronomers 
today in two significant ways. First, the signals are added instead of cross correlated (i.e., 
it is an additive interferometer) and, secondly, the two feeds involve different mix-down 
(i.e. Local Oscillator, or 'LO') frequencies. Therefore, when the signals are combined a 
beat signal (with frequency about 500 kHz) results. The end result, however, is that the 
response of the VSRT interferometer mimics, in many ways, a standard cross-correlation 
interferometer, as we show below. 

First consider the radiation from a single source entering the two feeds and the detected 
power that exits the square-law detector. We'll denote the baseline distance between the 
feeds as 'b,' and, assuming that the distance of the source is much larger than b, we assign 
the position of the source by its direction angle, 9, relative to the mid-plane between the 
feeds. The electric field entering each feed is 

E a = Eq cos(27rz/t) and = Eq cos(2-7rz/t — <fi), (5) 

where the phase difference, 0, is due to the extra path length to the second antenna and is 
given by 

b sin 9 

= 2vr^— . (6) 
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The signals are then mixed down with frequencies v a and z/&, added, and squared with only 
low-frequencies passing through. The end result of all these steps is^ 



V* = V 2 cos[2vr(z/ b - v a )t + 0] . (7) 

This, simply, is the beat signal. For simplicity of notation, we will write the amplitude of 
the beat signal as 'S', and so the last equation becomes 

S = S cos[27r(z/ fe — u a )t + </>]. (8) 

Now consider the total power in the beat signal when there are N sources. In general, the 
sources are incoherent and so we must add the powers. The total power in the beat signals 
of N sources, then, is given by 

St = ^ S k cos[27r(i/ 6 - v a )t + fc ], (9) 

k 

where Sk is the power in the beat signal from the kth source and <pk is the phase delay to 
the second feed from the kth source. This is mathematically identical tcP 



Vb — v a )t)- (10) 

Note that this, again, is the beat signal, but with an amplitude that is modified by the factor 
with the radical. Note also that the amplitude of the beat signal in Eq. [TU] is identical to 
Eq. |3j We see, therefore, that the power of the beat signal with the VSRT interferometer is 
identical to the amplitude of the complex visibility for that baseline for any arrangement of 
sources. 



IV. THE UNDERGRADUATE LABS 



Below we discuss a sequence of four labs designed to impart an intuitive sense of how data 
from an array of antennas can produce images of radio sources. Depending on the goals and 
conditions of the class the instructor could choose to skip the last lab and still successfully 
convey the basics of the relation between the source structure and the data between pairs 
of antennas. 
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FIG. 1: Undergraduate students from the Union College radio astronomy independent study group 
work with the VSRT interferometer. Two feeds taken from TV satellite dishes, are aimed at two 
CFL's, which act as radio sources. The feeds shown are "triple" DirecTV feeds for receiving 3 
satellites at 101, 110, and 119° West in geostationary orbit. The lowest feed of each triple as shown 
is the active feed. The nominal local oscillator frequency is 11.25 GHz for reception fo the 11.7 1 
12.2 GHz band. The feed polarization is left circular (LCP) which becomes right circular (RCP) 
upon reflection from the dish when used for TV. The distance between the feeds is determined 
by reading their positions relative to the meter stick below them. The output spectrum from the 
interferometer is displayed on the laptop screen, and the data files are recorded. The CFL's, which 
are easily movable, are located two-meters from feeds, which can be moved to assorted separation 
distances. 

For the sake of simplicity the lab set-ups are two-dimensional. Maps of the sources, then, 
are merely plots of intensity vs. position in the horizontal direction. In Figure [T] an example 
of the set-up is shown. 

Complete instructions for all of these labs is also available at 
https: / / wwwl.union.edu/marrj / radioastro /labfiles.html. 
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A. Analysis Software 



To facilitate the analysis of the data for these labs we have produced a package of java 
programs into which the output data from the VSRTI inteferometer control program are 
easily fed. These programs also have links to the lab instructions and can produce overlays 
of theoretical models with adjustable parameters. 

We have called this package "VSRTLPlotter" . A free and downloadable zip file of the VS- 
RTLPlotter package is located at https://wwwl.union.edu/marrj/radioastro/labfiles.html. 

B. Lab #1: The VSRT Primary Beam 

In the first lab, the students use the VSRT interferometer to observe a single CFL at 
varying angular positions and discover that the detected power depends on the position of 
the CFL. They learn about the "primary beam" of the individual antennas and that when 
sources are not at the center of the primary beam, the detected power is decreased by the 
primary beam factor. This is an exercise that's relevant to single-dish radio astronomy as 
well, since the beam size of the antenna determines the resolution of the telescope and the 
step size that a single-dish telescope must move to map a source. Additionally, the size of the 
primary beam places an effective upper limit on the maximum field of view in an aperture 
synthesis observation with a single pointing. Knowledge of the primary beam pattern is also 
important to the analysis in the fourth lab. 

The VSRTLPlotter program enables the students to overlay a theoretical beam plot and 
obtain a fitted value of the diameter of each antenna. 

In Figure[2]we show the beam pattern plot displayed in VSRTLPlotter with data obtained 
by undergraduate students in the radio astronomy class at Union College. 

C. Lab #2: The Visibility Function of a Single Resolved Source: 

Using a single CFL, the students are instructed to vary the separation distance between 
the feeds and plot the measured power vs. baseline distance and so are introduced to the 
"Visilibity function." Figure [3] displays data obtained by undergraduate students. The stu- 
dents see that the detected signal decreases with baseline distance. They then place two 
CFL's side-by-side, to simulate a source twice as large and discover that the decrease in the 
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FIG. 2: A plot of detected power vs. angular position of a CFL, revealing the VSRT's primary beam. 
The model overlay yields a measure of the effective diameter of the antenna feed of approximately 
3.8 cm. 

Visibility function is faster when the source is wider. Without delving into the mathemat- 
ics of the Visibiliity function, the students gain the intuition that a plot of interferometer 
response vs. baseline distance contains information about the angular size of the source. 

With the overlay model option in VSRTLPlotter, the students can fit the observed data 
and find a best-fit measure of the angular size of the source and see that it agrees with that 
expected. 

D. Lab #3: Visibility Function for a Pair of Sources 

Here, the students observe two CFLs separated by a given distance and again observe 
with the feeds at varying separations. The students discover that when observing a double 
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FIG. 3: The detected power when observing a single CFL vs. the separation distance between the 
feeds. The plot shows that with increasing baseline lengths more of the flux of the source gets 
resolved out. The smooth curve represents the model in which the brightness and angular width 
of the source is adjusted to fit. The same measurements can be made with two CFL's positioned 
side-by-side to model a single source twice as wide. Students then discover that the detected power 
decreases twice as fast. 

source the visibility function oscillates. Data for this lab obtained by students are displayed 
in Figure |4| By changing the separation distance of the CFLs and re-observing with the 
same set of baselines they discover that the oscillation length of the Visibility function is 
inversely related to the separation distance of the two sources. 
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FIG. 4: The detected power with a pair of CFLs vs. baseline length. The detected signal has an 
oscillating dependence on baseline length. The overall decrease with baseline is due to the single- 
resolved source pattern of Lab #2. The experiment can be repeated for different CFL separations 
to demonstrate that the periodicity of the oscillating function is inversely related to the separation 
angle of the sources. 

Again, with VSRTLPlotter, the students can fit a model to the data to obtain measures 
of the angular separation of the sources and compare with their set-ups in the lab. 
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E. Lab #4: Examination of Interferometer "Fringe Pattern" 

In this final lab, the students discover the "fringe pattern" of an interferometer. Although 
for some classes the principle contained in this lab may seem too esoteric, for students with 
strong backgrounds likely to continue on in astronomy, or for students starting research 
projects in radio astronomy, this lab exercise provides a visual demonstration of a funda- 
mentally important concept that newcomers to aperture synthesis commonly find confusing. 

When an interferometer is used to observe a celestial point source the direction of the 
source changes continuously as the Earth rotates causing the resulting detected power vs. 
time to oscillate quasi-sinusoidallyPEl The "fringe frequency" depends on the source position 
and is greatest when the source is near transit for an East- West baseline. One can also speak 
of the "fringe pattern," in which the detected signal is mapped as a function of the source's 
position in the sky. The fringe pattern is important to interferometery in three ways. First, 
fitting this function to the detected power vs. time yields an accurate measure of a source's 
position in the sky. Secondly, this function is used to calibrate interferometer baselines by 
observing a bright point source whose position is well known and fits the fringe function. 
And, thirdly, this function is used to determine when a faint source is being detected. If the 
power has a time dependence that fits the fringe function then the observer knows that it is 
due to a celestial source. 

To reproduce the fringe pattern with the VSRT interferometer, since the VSRT interfer- 
ometer cannot measure the fringe phases directly, one must place and fix one CFL at 9 = 
as a reference source, while the fringes due to a second CFL is detected. Additionally, since 
the source doesn't naturally move relative to the baseline, the students must manually move 
the second CFL to various positions. While leaving the feeds stationary at a given baseline 
and one CFL at the center position, the students record data with the second CFL located 
at a range of angles at a constant distance to the mid-point between the feeds. The detected 



power, then, is given by Eq. 10 with just two sources, 9i = and 62 = 9. 



S = y/Sf + SI + 25i5 2 cos(0) (11) 

A data set for this lab obtained by undergraduate students is shown in Figure [5] The 
amplitude of the oscillations decrease toward larger angles because the CFL moves out of 
the center of the primary beam. 
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FIG. 5: A plot of the detected power vs. angular position of a CFL, while another CFL stays fixed 
at 9 = 0. The resultant plot reveals the "fringe pattern,", although with an outer envelope due 
to the primary beam pattern. The model overlay, obtained for a baseline of 12.5 cm, more clearly 
shows how the fringe rate depends on the source position. Note that the distance between peaks 
when the source is near = is 0.2 radians, while at larger angles the peaks occur at 0.4 and 0.65 
radians. 

V. CONCLUDING REMARKS 

A group of undergraduate students undertaking an independent study on radio astronomy 
at Union College in Fall 2010 completed the set of labs discussed above. In an independent 
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study these students would not have succeeded in digging through the many pages of high 
level math needed to understand the basics of aperture synthesis and so these students would 
have completed their study considering only single-dish observing methods. Instead, these 
students ended the term with first-hand experience that demonstrated that obtaining data 
between pairs of antennas at many different separation distances leads to a function from 
which one could infer the size of a single source or the separation of a pair of sources. With 
a little further discussion, and possibly demonstrations using images from actual aperture 
synthesis data, an understanding of the concept that the structural details of even a complex 
source can be extracted from such data straightforwardly using established algorithms is 
obtained. 
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